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1.0  INTRODUCTION 


This  Annual  Scientific  Report  summarizes  research  accomplished  during 
the  second  year  (January  83  to  January  84)  of  AFOSR  Grant  No.  AFOSR-82- 
0041A.  The  objective  of  this  effort  is  to  determine  the  manner  in  which 
"small”  fatigue  cracks  initiate  at  notches,  extend  by  cyclic  loading, 
interact  with  adjacent  flaws,  and  coalesce  into  a  single  dominant  crack 
which  controls  final  fracture.  The  desired  product  is  a  predictive 
scheme  capable  of  analyzing  the  early  stages  of  fatigue  crack  growth 
which  are  characterized  by  the  growth  and  link-up  of  small  cracks. 

Research  toward  this  goal  is  directed  at  the  following  tasks. 

I.  Crack  Growth  Predictive  Algorithm.  A  computer  program  is  being 
developed  to  predict  the  growth  and  coalescence  of  multiple 
cracks  located  at  notches. 

11.  Crack  Interaction  Analysis.  Stress  intensity  factors  solutions 
are  computed  for  multiple  cracks  located  at  an  open  hole.  These 
solutions  are  required  for  the  multiple  crack  growth  algorithm. 

III.  Crack  Coalescence  Experiments.  Fatigue  tests  are  being  conducted 
with  multiply  cracked  specimens  to  provide  a  data  base  to  evaluate 
the  predictive  model.  Initially  the  model  will  be  verified  with 
"large  crack"  results  directed  toward  coalescence  aspects  of  the 
problem.  Subsequent  experiments  will  focus  on  coalescence  of 
"small"  cracks. 

IV.  Characterization  of  Small  Cracks.  This  phase  of  the  effort  is 
directed  toward  the  growth  and  coalescence  of  physically  small 
cracks.  These  data  will  be  used  to  characterize  the  growth  of 


“small"  flaws,  which  are  expected  to  behave  differently  than 
“large"  cracks. 

Additional  details  of  these  goals  and  progress  to  date  are  described  in 
the  remaining  sections  of  this  report. 


2.0  SUMMARY  OF  PROGRESS  AND  ACCOMPLISHMENTS 

2.1  Task  I  Progress  -  Crack  Growth  Algorithm 

A  computer  program  has  been  written  to  analyze  the  multiple  crack 
configurations  shown  in  Fig.  1.  Various  combinations  of  surface  and/or 
corner  flaws  are  assumed  located  along  the  bore  of  a  circular  hole.  It 
is  assumed  that  linear  elastic  fracture  mechanics  are  valid  and  that 
crack  growth  rates  are  controlled  by  the  stress  intensity  factor  K. 
Stress  intensity  factors  are  estimated  at  the  major  and  minor  crack 
axes  (locations  1  to  6  in  Fig.  1)  by  modifying  the  solutions  for  single 
cracks  In  a  large  body  [1]  with  an  interaction  factor  [2]  developed  here 
for  coalescing  cracks.  (The  interaction  solution  is  discussed  in  the 
Task  II  progress  section).  The  crack  tips  are  allowed  to  grow  indepen¬ 
dently  at  their  major  and  minor  axes  without  specifying  crack  shape 
changes.  The  multi-degree  of  freedom  program  predicts  the  growth  of 
the  separate  initial  cracks  and  their  coalescence  into  single  surface 
corner,  or  through-the-thickness  flaws  as  shown  in  Fig.  2.  Sample 
predictions  and  comparisons  with  test  data  are  described  in  the  Task 
III  progress  section. 

2.2  Task  11  Progress  -  Crack  Interaction  Analysis 

A  key  element  of  the  overall  predictive  scheme  is  determining  the 
interaction  between  adjacent  cracks  prior  to  their  actual  coalescence. 
The  presence  of  a  nearby  flaw,  for  example,  causes  in  increased  crack 
growth  rate  as  the  two  cracks  approach  each  other.  Although  a  few 
authors  have  studied  the  magnification  In  stress  intensity  factors 
due  to  coalescing  cracks  [3-6],  solutions  were  not  available  for 
multiple  cracks  at  notches.  Thus,  one  phase  of  the  current  effort  is 


directed  toward  obtaining  stress  intensity  factor  solutions  for  coalesc¬ 
ing  cracks  at  fastener  holes,  a  geometry  which  is  also  studied  experi¬ 
mentally  in  the  current  program. 

The  three-dimensional  finite  element-altering  method  [7-9]  was 
used  to  compute  stress  intensity  factors  for  symmetric  corner  cracks 
located  at  opposite  sides  of  a  plate  containing  a  hole  loaded  in  remote 
tension.  Although  the  analysis  was  limited  to  symmetric  corner  cracks, 
several  flaw  shapes  and  si;es  wore  considered.  By  comparing  the  double 
crack  results  with  corresponding  solutions  for  single  cracks,  it  was 
possible  to  determine  the  effect  of  the  second  crack  on  the  stress 
intensity  factor  as  the  two  flaws  grow  toward  each  other  and  eventually 
touch. 

Additional  details  of  the  crack  interaction  analysis  are  given  in 
Refs.  2  and  10.  Those  results  have  been  described  by  a  crack  inter¬ 
action  factor  which  has  been  incorporated  in  the  computer  program 
described  in  Task  1  to  account  for  crack  coalescence  in  the  overall 
prediction  scheme.  Experimental  verification  of  the  crack  coalescence 
predictions  are  described  in  the  following  section. 

2.3  Task  Ill  Progress  -  Crack  Coalescence  Experimental  Evaluation 

This  phase  of  the  effort  is  aimed  at  evaluating  the  numerical 
scheme  described  in  Task  I  for  predicting  the  initial  growth  and 
coalescence  of  multiple  cracks  located  at  notches.  In  order  to 
separate  crack  coalescence  aspects  from  "small  crack"  phenomena,  the 
initial  work  involved  experiments  with  coalescence  of  relatively 
"large"  cracks.  Current  work  with  the  small  crack  problem  is  described 
separately  in  the  Task  IV  section. 


As  suii«narized  in  References  11  and  12,  verification  of  the  predic¬ 
tive  crack  growth  model  has  included  comparison  with  several  sets  of 
experimental  data.  Both  single  cracks  at  holes  and  multiple  coalescing 
flaws  have  been  studied.  The  multiple  crack  configurations  include 
coalescing  cracks  at  open  holes  in  polymethylmethacylate  (PMMA)  plates 
(transparent  polymer  specimens  which  allow  direct  observation  of  the 
crack  plane),  PMMA  bend  specimens,  and  metal  edge  notch  specimens. 

Excellent  results  [11,12]  have  been  obtained  for  large  corner 
cracks  which  coalesce  along  the  bore  of  an  open  hole  in  a  PMMA  tension 
specimen  (figure  lc  configuration)  and  for  small  surface  cracks  which 
occur  along  a  semicircular  edge  notch  in  metal  specimens.  The  PMMA 
hole  crack  data  were  generated  during  the  present  program  [13],  while 
the  metal  data  were  obtained  from  another  investigation  [14]. 

The  metal  results  described  in  an  earlier  progress  report  [12] 
were  obtained  by  Or.  A.  Thakker  with  Wa spa  Hoy  specimens  shown  in 
Figure  3.  Additional  test  results  have  been  obtained  by  Thakker  for 
titanium  specimens  (Ti-6-2-4-6),  and  crack  growth  predictions  were 
performed  here.  Typical  predictions  for  the  titanium  small  crack 
coalescence  tests  are  given  in  Figures  4-7,  while  Fig.  8  summarizes 
predicted  and  test  lives  for  the  titanium  and  Waspalloy  experiments. 

The  tests  were  conducted  with  the  edge  notched  specimen  geometry  shown 
in  Figure  3.  Crack  lengths  were  measured  from  fracture  surface 
markings  made  by  heat  tinting  techniques.  (The  specimens  were 
periodically  heated  to  oxidize  the  crack  surfaces).  Note  that  the 
predictive  model  generally  gives  an  excellent  estimate  for  total 
specimen  life. 
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2 . 4  Task  IV  Progress.  -  Characterization  of  Small  Cracks 

It  is  commonly  known  that  physically  small  cracks  may  grow  faster 
than  larger  flaws  subjected  to  similar  stress  intensity  factor  levels. 

The  assumptions  regarding  continuum  behavior  and  linear  elastic  material 
response  lose  validity  as  crack  lengths  become  increasingly  small.  The 
problem  is  further  complicated  by  the  presence  of  adjacent  flaws  which 
often  initiate  simul taneously  and  which  interact  prior  to  the  coalescence 
process.  The  preceeding  sections  have  described  tasks  directed  toward 
the  coalescence  problem.  This  prior  work  has  concentrated  on  relatively 
large  flaws  in  order  to  separate  small  crack  behavior  from  coalescence 
phenomena.  This  section  briefly  outlines  work  aimed  at  the  "small  crack 
problem." 

A  series  of  fatigue  tests  have  been  conducted  with  V-notched  PMMA 
bend  spe  mens.  Natural  fatigue  cracks  were  allowed  to  develop  along 
the  V-notch  and  to  coalesce  into  a  single  flaw.  Time  lapse  photography 
was  used  to  record  and  measure  the  individual  and  coalescing  flaws. 

Flaw  sizes  (and  shapes)  less  than  0.005  Inches  are  routinely  observed 
and  measured.  Typical  results  are  given  in  References  12  and  13. 

In  a  few  tests,  the  multiple  cracks  quickly  joined  to  form  short 
through-the-thlckness  flaws.  Stress  intensity  factors  have  been 
estimated  for  these  small  cracks,  and  crack  growth  rates  are  compared 
with  baseline  data  for  larger  through -era eked  specimens  in  Fig.  9. 

These  preliminary  results  indicate  that  small  cracks  grow  faster  than 
large  flaws  (at  similar  sK  levels),  as  typically  observed  in  structural 
metals. 

Analysis  of  these  tests  Is  complicated  by  the  stress  field  present 
at  the  sharp  V-notch.  Additional  analysis  of  the  notch  specimen  are 
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being  conducted  to  verify  the  snail  crack  behavior  shown  in  Fig.  9,  and 
to  allow  predictions  for  the  growth  and  coalescence  of  multiple  surface 
cracks  at  notches.  It  is  recognized  that  this  goal  could  be  complicated 
by  nonlinear  plastic  zone  effects  which  may  limit  the  LEFM  approach  for 
*he  snail  crack  tests. 
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3.0  PLANNED  RESEARCH 

This  section  briefly  summarizes  research  planned  for  the  remainder 
of  the  proposed  effort.  For  convenience,  discussion  is  organized  around 
the  four  tasks  described  in  the  previous  sections. 

3 . 1  Task  i  pi  a jn>_  Crack  Growth  Algorithm 

The  current  crack  growth  model  predicts  growth  and  coalescence  of 
multiple  cracks  located  at  circular  holes.  Linear  elastic  fracture 
mechanics  are  assumed  to  be  valid  as  stress  intensity  factors  are 
used  to  compute  crack  growth  rates. 

The  multi-degree  of  freedom  predictive  program  has  been  written 
in  a  modular  manner,  so  that  modifications  may  be  made  as  required. 

It  is  expected  that  new  stress  intensity  factor  solutions  will  be 
added  to  analyze  alternate  specimen  geometries  (e.g.  surface  cracked 
bend  and  semi-circular  notched  bend  specimens).  In  addition,  should 
subsequent  work  indicate  that  LEFM  does  not  apply  to  the  small  crack 
tests,  an  alternate  crack  growth  parameter  (to  AK)  may  be  incorporated 
in  the  predictive  scheme.  Current  literature  on  nonlinear  crack  growth 
(e.g.  Ref.  15)  will  be  studied  as  appropriate. 

3 . 2  Task  II  PI a ns  -  Crack  Interaction  Analysis 

Currently,  the  Newman/Raju  [1]  stress  intensity  factor  solutions 
for  cracked  holes,  modified  with  a  crack  Interaction  factor  [2]  to 
account  for  multiple  cracks,  form  the  basis  for  the  original  stress 
intensity  factor  solutions.  Although  these  stress  intensity  factor 
solutions  were  derived  for  cracked  holes  in  large  plates,  they  have 
been  successfully  modified  with  a  "correction  factor"  [11]  and  applied 
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to  semicircular  notches  loaded  in  remote  tension  (as  required  for  the 
crack  growth  predictions  shown  in  Figs.  4-7). 

It  will  be  necessary  to  estimate  stress  intensity  solutions  for 
the  other  specimen  geometries  employed  in  the  experimental  effort.  These 
other  specimen  types  include  the  semicircular  edge  notch,  V-notch,  and 
keyhole  specimens  shown  In  Fig.  10.  In  addition,  a  surface  cracked  bend 
specimen  is  also  being  investigated.  All  of  these  specimens  are  loaded 
in  four-point  bending  In  a  manner  which  provides  direct  observation  of 
the  crack  plane.  Although  these  specimens  are  readily  suited  for  the 
experimental  work,  K  solutions  are  not  available  for  the  notched  geo¬ 
metries,  although  Newman  and  Raju  provide  results  for  the  surface  crack 
loaded  in  bending  [17].  Stress  Intensity  factors  will  be  estimated  for 
the  notched  specimens  either  by  the  "correction  factor"  approach  success¬ 
fully  employed  for  the  semicircular  edge  notch  tension  specimen  [11]  or 
by  employing  the  "crack  face  pressure"  stress  Intensity  solutions  given 
in  Reference  16.  The  semicircular  notch  bend  specimen  may  be  examined 
in  more  detail  with  the  finite  element  alternating  method  employed  in 
Reference  2.  The  goal  here  Is  to  obtain  reasonable  stress  Intensity 
factor  estimates  for  the  experimental  effort,  and  not  to  conduct  an 
in-depth  stress  Intensity  factor  analysis  program  for  various  crack 
configurations. 

Currently  the  interaction  factors  developed  for  symmetric  corner 
cracks  at  holes  [2,10]  are  used  to  predict  the  increased  growth  rates 
as  two  cracks  approach  each  other.  The  present  analysis  assumes  the 
cracks  are  symmetric,  located  at  open  holes,  and  lie  in  the  same 
plane.  The  results  have  worked  well  for  cracked  holes  and  semicircular 
notches  loaded  in  tension.  As  experimental  results  are  obtained  for 
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nonsymmetric  and  slightly  out-of-plane  flaws,  it  may  be  necessary  to 
modify  the  interaction  solutions.  Again,  potential  modifications  will 
be  "engineering  type"  estimates,  since  a  rigorous  three-dimensional 
analysis  of  nonsymmetric,  out-of-plane  cracks  is  beyond  the  scope  of 
the  present  investigation. 


3.3  Task  III  Plans  -  Crack  Coalescence  Experimental  Evaluation 


The  goal  here  is  to  further  evaluate  the  multiple  crack  growth  and 
coalescence  scheme  through  comparison  with  additional  experimental  data. 
To  date,  successful  predictions  have  been  obtained  for  symmetric 
coalescing  corner  cracks  at  holes  In  PMMA  plates  [11,12],  and  for  small 
surface  flaws  at  semicircular  edge  notches  In  Waspalloy  [12]  and 
titanium  tension  specimens  (Figs.  3-8).  In  addition,  experimental 
results  [13]  have  been  obtained  for  coalescing  flaws  in  PMMA  bend 
specimens  (both  the  semicircular  and  V-notched  configurations  shown 
In  Fig.  10).  As  mentioned  In  Section  3.2,  efforts  are  underway  to 
estimate  stress  Intensity  factors  for  these  two  notched  bend  specimens. 
When  completed,  those  solutions  will  be  incorporated  In  the  life 
prediction  program  and  predictions  made  for  these  PMMA  test  results. 

To  expand  the  evaluation  to  a  larger  data  base,  tests  are  being 
conducted  with  polycarbonate  specimens.  Polycarbonate  is  also  a 
transparent  polymer,  but  exhibits  considerably  more  ductility  than 
PMMA.  Baseline  fatigue  crack  growth  data  have  been  collected  for  the 
polycarbonate  material  and  are  given  In  Fig.  11.  Note  that  stable 


fatigue  crack  growth  rates  (da/dN)  are  obtained  at  aK  levels  of  3,000 


1/2 

psl-in  '  for  the  polycarbonate,  while  as  shown  In  Fig. 
PMMA  Is  on  the  order  of  1000  psl-in1^.  (Specific  K, 


9,  K.  for  the 


measurements  will 


will  be  performed  for  both  the  PMMA  and  polycarbonate  test  material, 
along  with  stress/strain  curves. 

Polycarbonate  crack  coalescence  tests  are  being  conducted  with 
cracked  hole  tension  specimens  and  with  a  surface  cracked  bend  geometry. 
Experiments  have  also  been  conducted  with  the  notched  bend  specimens 
shown  in  Fig.  10.  The  cracked  hole  tension  specimen  shown  in  fig.  12 
has  been  analyzed  by  several  investigators  [1,2,7,16],  and  K  solutions 
are  well  established  for  many  flaw  shapes.  In  addition,  existing  finite 
element-alternating  method  computer  codes  [2,16]  may  be  applied  to 
analyze  new  flaw  shapes  of  interest  to  the  current  effort. 

The  flawed  hole  specimen  is  somewhat  awkward  for  experimental  work, 
however,  since  grips  must  be  glued  to  the  specimen  ends  and  the  crack 
plane  viewed  with  a  mirror.  More  significant,  however,  is  the  fact 
that  relatively  large  hole  sizes  are  employed.  Although  this  is  not 
a  disadvantage  for  coalescence  tests  with  "large"  cracks,  where  large 
sta.ter  flaws  can  be  employed,  the  specimen  is  not  as  suited  to  small 
crack  studies.  Since  the  hole  radius  is  farily  large,  natural  flaws 
would  not  be  expected  to  initiate  in  the  same  plane,  and  the  coales¬ 
cence  study  would  be  complicated  by  out  of  plane  cracks,  (see  Figure 
13)  Although  it  will  eventually  be  necessary  to  study  the  influence 
of  out-of-planeness  on  crack  coalescence,  it  appears  preferable  to 
initially  study  flaws  which  grow  and  coalesce  in  the  same  (or  nearly 
same)  plane. 

In  order  to  employ  a  notch  geometry  with  a  smaller  radius  of 
curvature,  and,  thus,  confine  crack  initiation  and  growth  to  the 
same  plane,  the  V-notch  and  keyhole  specimens  shown  in  Fig.  10  were 
examined.  Both  specimens  were  subjected  to  four-point  bending  to 


simplify  load  application  and  viewing  of  the  crack  plane.  All  of 
these  specimens  have  worked  well  for  experimental  studies  on  the 
initiation  of  small  fatigue  cracks.  Unfortunately,  however,  the  notch 
tip  stresses  significantly  complicate  analysis  of  the  experimental 
results. 

Current  work  has  focused  on  an  unnotched  bend  specimen.  Small 
starter  slots  are  placed  adjacent  to  each  other  in  the  desired  location 
on  the  tension  surface  of  a  beam  with  a  rectangular  cross  section.  The 
specimen  is  annealed  to  remove  residual  stresses  introduced  by  the 
preflaws,  and  the  specimens  are  cycled  in  four-point  bending.  Small 
fatigue  cracks  develop  around  the  preflaws,  and  grow,  eventually  coal¬ 
escing  into  a  single  surface  or  corner  crack.  The  Newman  and  Raju 
stress  intensity  factor  solution  for  single  surface  cracks  loaded  in 
bending  [17]  will  be  modified  with  the  crack  interaction  factor  as 
before,  and  incorporated  in  the  life  prediction  scheme  to  predict 
specimen  life. 

In  summary,  then,  work  on  this  task  will  concentrate  on  conducting 
more  crack  coalescence  tests  with  polycarbonate  and,  perhaps,  PMMA 
specimens.  Multiply  cracked  hole  specimens  will  yield  data  on  coal¬ 
escence  of  fairly  large  flaws.  It  is  planned  to  look  at  nonsymmetric 
corner  cracks,  and  surface/corner  crack  combinations  which  have  not 
been  studied  to  date.  This  specimen  closely  matches  the  boundary 
conditions  assumed  in  the  current  analysis  scheme  (Task  I),  and  is  quite 
suited  for  crack  coalescence  studies.  Initiation  and  coalescence 
studies  with  smaller  fatigue  cracks  will  employ  the  bend  specimen. 

Both  polycarbonate,  and,  perhaps,  few  PMMA  tests  will  be  conducted. 
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As  mentioned  previously,  this  specimen  yields  small,  natural  cracks, 
which  lie  in  nearly  the  same  plane. 

Although  the  original  proposal  suggested  a  small  pilot  study  with 
metal  specimens,  it  is  not  believed  that  these  tests  are  necessary  now, 
in  view  of  the  excellent  results  which  have  been  obtained  with  the 
Waspalloy  (Figs.  18  to  23  in  Ref.  12)  and  Titanium  (Figs.  4  to  8  in 
this  document)  specimens  tested  by  Thakker  [14].  It  is  felt  more 
appropriate  to  concentrate  on  further  tests  with  PMMA  and  polycarbonate 
specimens,  where  crack  shapes  and  sizes  are  more  accurately  measured, 
rather  than  diluting  the  effort  on  a  small  scale  program  with  metal 
specimens. 

3.4  Task  IV  Plans  -  Characterization  of  Small  Cracks 

The  goal  here  is  to  characterize  the  growth  of  very  small  fatigue 
crack  sizes  to  determine  if  they  grow  differently  than  large  flaws. 
Efforts  will  focus  on  the  growth  of  single  cracks  (or  widely  separated 
multiple  flaws)  so  that  small  crack  behavior  can  be  separated  from 
coalescence  phenomena.  Since  the  cracks  must  have  small  physical 
dimensions,  they  will  be  initiated  "naturally"  at  local  stress  con¬ 
centrations.  It  is  important  that  the  starter  notches  have  well 
defined  dimensions,  so  that  the  local  stress  level  can  be  determined 
and  appropriate  crack  parameters  computed.  Initially,  for  example, 
crack  growth  rates  will  be  correlated  with  the  cyclic  stress  intensity 
factor. 

Currently  the  crack  photographs  are  analyzed  by  projecting  the 
35  mm  negatives  onto  a  screen  and  measuring  the  enlarged  image  (approxi¬ 
mately  30  x)  to  the  nearest  mm.  Since  this  method  is  somewhat  crude 


and  cumbersome,  a  film  positioner  has  been  built  to  measure  the  film 
strips  under  a  microscope.  The  film  holder  consists  of  two  miciometer 
positioners  mounted  perpendicular  to  each  other.  The  micrometer  heads 
are  calibrated  in  0.0001  in.  increments,  allowing  x-y  coordinates  to 
be  located  by  cross  hairs  in  the  microscope  and  measured  quite  accurately 
on  the  film  strip.  Steps  have  been  taken  to  further  automate  the  film 
reader  with  electrical  readout  to  display  crack  lengths  automatically 
in  digital  form,  and,  perhaps,  be  recorded  directly  for  computer  analysis. 
It  is  felt  the  new  film  reader  will  significantly  improve  accuracy  and 
reduce  the  effort  required  to  measure  the  crack  lengths. 

In  summary,  this  task  will  focus  on  measuring  crack  growth  rates 
for  small  flaw  sizes  in  both  polycarbonate  and  PMMA  specimens.  Stress 
intensity  factors  will  continue  to  be  used  to  characterize  crack 
behavior,  unless  the  results  indicate  LEFM  is  no  longer  valid.  In 
that  case,  an  alternate  crack  tip  parameter,  such  as  a  strain  intensity 
factor  or  the  J  integral,  may  be  examined.  If  necessary,  the  multicrack 
growth  and  coalescence  prediction  scheme  will  then  be  modified  to 
incorporate  the  new  crack  growth  characterization  term. 
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Figure  1  Schematic  of  multiple  cracks  located  along  bore  of  a  hole 
in  a  large  plate.  Hole  diameter  is  D,  plate  thickness  is  t 
and  remote  tensile  stress  is  applied  in  z  direction. 
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Figure  3  Drawing  of  test  specimen  showing  specimen  dimensions 
and  location  of  two  semicircular  edge  notches. 


arison  of  predicted  and  experimental  growth  and  coalescence  of 
surface  cracks  at  semicircular  notch  In  a  titanium  specimen 


TEST  10B  :  EDGE  NOTCHED  TENSILE  COUPON 
TWO  SURFRCE  CRACKS 


TEST  12B  :  EDGE  NOTCHED  TENSILE  COUPON 
TWO  SURFACE  CRACKS 


arison  of  predicted  and  experimental  growth  and  coalescence  of 
surface  cracks  at  a  semicircular  notch  in  a  titanium  specimen 


Comparison  of  predicted  and  experimental  growth  and  coalescence  of 
two  surface  cracks  at  a  semicircular  notch  in  a  titanium  specimen 
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PREDICTED  VS  TESTED  CYCLES  TO  FRILURE 


Figure  8  Comparison  of  total  crack  growth  life  predicted  by  current 
model  with  results  of  crack  coalescence  experiments 
conducted  by  Thakker  with  Waspalloy  and  Titanium  test  specimens. 
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Figure  10  Schematic  drawing  of  semiciruclar,  V-notch,  and  keyhole 
specimens  showing  4-point  bend  loading  and  placement  of 
viewing  camera  (not  to  scale). 
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Figure  11  Baseline  fatigue  crack  growth  rate  da/dN  versus  AK 
for  "large"  crack  polycarbonate  test  specimens 
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STRESS  INTENSITY  FACTORS  FOR  COALESCING 
AND  SINGLE  CORNER  FLAWS  ALONG 
A  HOLE  BORE  IN  A  PLATE 

B  J.  HEATH  and  A.  F.  GRANDT.  Jr. 

School  of  Aeronautics  and  Astronautics.  Purdue  University.  W.  Lafayette.  IN  47907,  U.S.A. 

Abstract— The  objective  of  this  paper  is  to  describe  the  effects  of  crack  interaction  on  stress  intensity 
factors  for  two  symmetric  coplanar  comer  flaws  located  along  a  hole  bore.  This  numerical  analysis 
employes  the  Finite  Element-Alternating  Method  to  determine  Mode  1  stiess  intensity  factors  for  single 
and  coalescing  corner  flaws.  Using  single  flaw  stress  intensity  factors  as  a  reference,  analysis  of  crack  size 
and  shape  effects  on  Jij  for  coalescing  corner  flaws  indicates  the  stress  intensity  factor  for  crack  points 
along  the  hole  bore  increases  as  the  crack  tip  separation  distance  decreases.  Interaction  effects  are  not 
enperienced  by  hole  bore  cr.ick  points  w  hen  the  crack  lip  separation  distance  is  equal  to  or  greater  than  half  of 
the  largest  corner  flaw  dimension. 


INTRODUCTION 

Naturally  occurring  fatigue  cracks  frequently  initiate  at  several  independent  points  along  the  bore  of 
fastener  holes.  These  individual  flaws  then  grow  and  coalesce  into  a  single  dominant  crack  which 
controls  final  fracture.  Prediction  of  this  coalescence  phase  of  crack  growth  requires  stress  intensity 
factor  solutions  for  the  individual  cracks  prior  to  their  “link  up"  into  a  single  flaw. 

It  is  known  that  the  stress  intensity  factor  for  a  crack  tip  is  influenced  by  the  close  proximity  of  an 
adjacent  crack  tip  [  1— 4J.  Crack  interaction  effects  on  Mode  I  stress  intensity  factors  (K,)  have  been 
studied  by  several  authors  11-4].  Various  2-dimensional  Kt  solutions  exist  for  interaction  of  coplanar 
through-the-thickness  cracksl3,4].  Kamei  and  Yokoburi[3],  present  a  2-dimensional  stress  intensity 
factor  solution  for  two  asymmetric  through-the-thickness  cracks  in  an  infinite  sheet,  while  Benthem  and 
Koiter  [4]  give  a  2-dimensional  K,  solution  for  two  symmetric  edge  cracks.  Murakami  and  Nemat- 
Nassar  [1]  present  an  approximate  3-dimensional  K,  solution,  obtained  by  the  Body  Force  Method,  for 
coplanar  surface  flaws  in  an  unbounded  solid.  None  of  these  solutions,  however,  consider  coalescence  of 
coplanar  flaws  located  along  the  bore  of  a  hole  in  a  finite  thickness  plate.  Since  multiple  fastener  hole 
cracks  often  occur  in  service  15,6],  the  present  paper  examines  crack  coalescence  along  a  hole  bore. 

The  objective  of  this  paper  is  to  describe  stress  intensity  factor  results  for  two  adjacent  cracks 
located  along  the  bore  of  a  hole  in  a  finite  thickness  plate.  Figure  1  shows  typical  plate  cross  sections  for 
the  symmetric  and  single  crack  configurations  considered.  Here  a  is  the  crack  dimension  measured  along 
the  hole  bore;  c  is  measured  along  the  free  surface;  T  is  the  plate  thickness  and  D  is  the  hole  diameter. 
The  plate  is  loaded  with  a  remote  stress  a  applied  perpendicular  to  the  crack  plane.  The  Finite 
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Fig.  I.  Schematic  drawing  of  crack  plane  showing  location  of  single  and  symmetric  comer  cracks  located  at 

bote  of  hole. 
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l  It  ;,it  n!  Alter rutlii..:  Method  (I  I  AM)  (7-  13]  is  us.  J  to  compute  Kt  for  s;.  n  metric  (coalescing)  and  single 
coiner  cracks  alorg  the  hole  bore. 

Results  are  given  in  both  tabular  and  graphical  form  for  crack  sizes  in  tht  range  0.25  <  a/7  <0.5  for 
crack  shapes  ale  =  1.11,  1.5,  2.0  and  3.0.  (As  described  in  earlier  work  17),  numerical  programming 
difficulties  present  analysis  of  quarter  circular  cracks,  ale  =  1.0,  with  the  present  computer  codes.) 
Crack  shape  and  si/e  effects  on  Kt  variation  around  the  crack  perimeter  are  presented  for  coalescing 
corner  flaws.  Two  and  3  dimensional  sti^ss  intensity  factor  solutions  for  other  coalescing  crack 
configurations  are  compared  with  the  stress  intensity  factors  generated  for  the  corner  cracked  hole 
geometry. 

NUMERICAL  APPROACH 

Stress  intensity  factors  for  coatescencing  cracks  along  a  hole  bore  arc  obtained  by  the  Finite 
Element- Alternating  Method  (FEAM).  The  present  parametric  study  makes  use  of  computer  codes 
developed  by  Smith  and  Kullgren [7].  The  original  codes  were  modified  here  to  treat  the  coalescing  crack 
problem. 

The  FEAM  calculates  Kt  based  on  approximate  surface  crack  boundary  conditions.  An  iterative 
superposition  of  two  solutions,  one  a  3-D  finite  element  solution  for  an  uncracked  body  under  prescribed 
surface  loadings,  and  the  second  a  stress  solution  for  a  flat  elliptical  crack  in  an  infinite  body  with 
nonuniform  prescribed  surface  pressure,  produce  these  approximate  surface  flaw  boundary  conditions. 
Earlier  applications  of  the  FEAM  to  various  other  cracked  hole  problems  are  described  in  Refs  [8-12], 

The  crack  coalescence  problem  studied  here  requires  modification  of  the  finite  element  mesh  and  the 
elliptical  crack  stress  solution  of  [7]  to  produce  the  symmetric  crack  finite  element  model  used  in 
conjunction  with  the  coalescence  option  of  the  FEAM  computer  code. 

Single  crack  model 

The  single  crack  model  is  a  semi-circular  plate  with  a  hole.  As  in  Fig.  2,  this  model  has  a  hole  radius 
R.  In  order  to  minimize  boundary  effects  in  the  region  of  the  crack,  the  same  outer  plate  radius  of  12  R 
(R  =  hole  radius)  used  by  Smith  and  Kullgren  [7]  is  used.  The  single  crack  model  employs  112  20-Node 
isoparametric  finite  elements.  The  hole  diameter  to  thickness  ratio  (D/T)  is  0.5.  The  modulus  of  rigidity 
for  the  plate  is  G  =  12  x  106psi  with  a  Poisson’s  ratio  of  v  =  0.3. 

Symmetric  crack  model 

Symmetric  interacting  cracks  were  studied  by  imposing  a  line  of  symmetry  along  the  x  =  0.5  face  of 
the  original  FEAM  model  with  D/T  =  1.0  as  shown  in  Fig.  2.  The  hole  radius,  outer  radius,  and  element 
number  are  equivalent  to  the  single  crack  model.  In  order  to  account  for  the  line  of  symmetry,  the 
symmetric  crack  model  has  a  D/T  =  1.0  rather  than  D/T  =  0.5  for  the  single  crack.  As  shown  in  Fig.  3, 
the  plane  of  symmetry  imposed  along  the  front  surface  of  the  plate,  however,  effectively  doubles  the 
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symmetric  crack  model  plate  thickness.  Thus,  the  symmetric  crack  model  D/T  equals  that  of  the  single 
crack  model  (D/T  -  0.5).  Figure  3  shows  how  an  infinite  plate  with  two  symmetric  cracks  along  the  hole 
bore  is  simulated  when  a  corner  flaw  is  applied  to  the  FEAM  symmetric  crack  model. 

DISCI  SSION  OF  RESl'LTS 

Stress  intensity  factors  were  generated  for  crack  shapes  ale  =  1.11,  1.5,  2.0  and  3.0  by  the  FEAM 
analysis.  F.ach  crack  varied  in  size  from  0.25  sa/T  <0.5.  The  stress  intensity  factors  for  the  corner 
flaw  configurations  considered  are  listed  in  Tables  1  and  2.  Table  1  lists  the  stress  intensity  factor 
variation  for  the  single  crack  in  15  degree  increments  of  the  elliptic  angle  defined  in  Fig.  1.  Stress 
intensity  factor  results  for  the  symmetric  crack  configuration  are  listed  in  Table  2. 

These  results  provide  a  parametric  study  of  size  and  shape  effects  on  the  stress  intensity  factor 
variation  around  the  crack  perimeter.  For  a  single  corner  crack  with  ale  -  1.11,  Fig.  4  shows  the 
variation  of  the  dimensionless  stress  intensity  factor  (K,la^D)  vs  normalized  elliptic  angle  as 

the  crack  size  (a/T)  increases.  Note  that  90  degrees.  Figure  5  shows  the  variation  of  K,la\/D 

over  the  same  a/T  range  for  two  symmetric  cracks  with  the  same  ale  =  1.11. 

Recalling  that  4>  -■  90'  cm  responds  to  the  hole  bore  crack  location,  notice  the  increased  magnitude  of 


Table  I  Dimensionless  stress  intensity  factors  foi  single  corner  crack  al  open  bolt  in  a  large  plate  loaded  in  remote 
tension  Results  gisen  as  a  function  of  crack  shape  rt'r.  crack  si/e  a>7.  and  position  alonp  crack  perimeter  defined  by 
parametric  angle  d  li  =  0'  al  fionl  surface  and  &  90*  al  holt  bore) 
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Fig.  5.  Dimensionless  stress  intensit)  factors  for  symmetric  double  cracl.  problem  as  function  of  dimen¬ 
sionless  parametric  angle  616 mu  (d  -  dmu  =  W  at  hole  bore)  and  crack  sire  alT  for  flaw  shape  ale  «  Ml. 
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fig  6  Comparison  of  dimensionless  stress  intensity  factors  for  sirg'e  and  symmetric  cracks  »i-F  aspect 
ratio  a'r  =  I  II  sho»  irg  nagr.if.calion  in  stress  irte-sit)  factor  at  hole  bort  location  16  90'l  as  cr^.l  sire 

alT  increases 
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Kilos/! )  along  the  hole  bore  for  the  symmetric  configuration  compared  to  the  same  location  for  the 
single  fl;iss  as  the  crack  si/e  is  increased.  Using  the  single  flaw  KJas/D  as  a  reference,  the  influence  of 
at T  on  Kifrrv'O  for  particular  fl.»ss  points  on  the  symmetric  model  crack  perimeter  is  shov.n  in  Fig-.. 
6-S.  f  or  small  crack  si/es  (0.2^  t-  alT  ^0.3)  shown  in  Figs.  6-  8,  the  single  and  coalescing  configurations 
differ  by  a  maximum  of  4T.  How  ever.  Figs.  6  and  7  show  a  significant  increase  in  the  symmetric  crack 
Kilos/ D  as  compared  to  the  single  flaw  K/os/D  when  the  crack  size  approaches  coalescence 
(a/T  =•  0.5). 

Folio's ing  Murakami  11),  a  crack  interaction  factor  y  is  defined  as  the  ratio  of  the  stress  intensity 
factor  for  two  cracks  to  the  stress  intensity  factor  for  a  single  crack.  Kamci  and  Yokoburi  [3J  have  given 
2-dimensional  stress  intensity  factor  solutions  for  two  asymmetrical  through-the-thickness  cracks  in  an 
infinite  elastic  sheet.  By  equating  the  crack  lengths  as  shown  in  Fig.  9(a),  and  normalizing  the  crack  tip 
stress  intensity  factor  expression  by  os/va,  the  interaction  factor  y  for  two  symmetric  through-the- 
thickness  cracks  is  obtained.  The  Bcnthem  and  Koiter  [4]  2-dimensional  stress  intensity  factor  solution 
for  two  symmetric  edge  cracks  (see  Fig.  9b)  may  be  used  to  compute  the  edge  crack  interaction  factor 
shown  in  Fig.  10.  Here,  the  double  edge-crack  solution  was  divided  by  the  stress  intensity  factor  given 
by  Harris  [14]  for  a  single  edge-cracked  sheet  which  doesn’t  bend. 
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Fig.  9.  Interactinp  crack  configurations  represented  by:  (a)  Kamei  (throuph-the-thiekness  center  cracks  in  a 
large  sheet),  (b)  Koiter  and  Bcr.them  (throuph  the-thickness  edge  cracks),  (c)  Murakami  and  Nemet-Nasser 
(surface  cracks  in  a  semi-infinite  solid),  and  (d)  present  study  (symmetric  corner  cracks  located  along  hole 

bore). 
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Fig.  10,  Comparison  of  effect  of  crack  spacing  f„p)a  on  stress  intensity  factor  for  various  2-  and 
3-dimensiona!  crack  configurations.  The  interaction  factor  y  is  defined  as  ratio  of  double  to  single  crack 

stress  intensity  factors. 


Three-dimensional  crack  point  interaction  factors  for  symmetric  semi-circular  surface  flaws  (Fig.  9c) 
arc  presenlcd  by  Murakami  and  Ncmat-Nassar  in  Ref.  [1],  These  coplanar  flaws  exist  in  an  unbounded 
solid.  The  present  FEAM  results  were  used  to  compute  3-dimcnsional  interaction  factors  for  symmetric 
corner  cracks  along  a  hole  bore  (Fig.  9d).  Figure  10  compares  these  2-  and  3-dimcnsional  crack  tip  interaction 
factors  as  a  function  of  dimensionless  crack  separation  distance  (t,rr/a).  The  through-the-lhickness 
interaction  factor  of  Kamei  (3)  represents  the  limiting  value  of  two  interacting  surface  cracks,  while  the 
Bcnfhcm  and  Koiter  |4]  edge-crack  factor  serves  as  a  limiting  value  for  the  present  study  of  interacting 
coiner  cracks.  Due  to  a  numerical  instability  in  the  FEAM  code,  values  of  crack  interaction  factors  at 
coalescence  (t„rfa  -  0)  are  not  obtainable. 

Notice  the  close  proximity  of  the  present  FEAM  results  for  corner  cracked  holes  to  the  Murakami  |1) 
interaction  factor  for  surface  cracks.  As  the  crack  separation  tirp  exceeds  approximately  half  of  the 
crack  length  dimension  a,  ( tupla  a  0.5),  the  FEAM  crack  interaction  factor  stabilize  at  unity,  indicating 
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that  the  crack  tip  located  along  the  hole  bore  is  no  longer  influenced  by  the  adjacent  crack.  For  a  given 
luJa,  the  maximum  difference  in  all  hole  bore  interaction  factors  presented  in  Fig.  10  is  10.7%.  Figure  10 
indicates  that  crack  size  (n/T)  dominates  the  increase  in  K(  due  to  crack  interaction  at  the  hole  bore, 
while  crack  shape  (ale)  has  only  a  slight  influence.  The  effect  of  crack  shape  on  the  hole  bore  interaction 
factors  is  shown  on  an  expanded  scale  in  Fig.  11. 

Figure  12  presents  the  interaction  factor  at  the  plate  surface  ( <j>  -  0).  The  factors  decrease  slightly 
below  unity  in  the  range  of  0  <  (/wp/fl)  s 0.5.  Figure  12  further  indicates  that  the  crack  is  free  of  influence 
from  the  adjacent  crack  when  the  hole  bore  separation  distance  is  greater  than  half  the  crack  length  a. 
As  the  separation  distance  decreases,  the  decrease  in  K(  at  the  plate  surface  is  considerably  less  than  the 
increase  in  Kt  at  the  hole  bore  location  (see  Fig.  10).  Thus,  the  effect  of  crack  coalescence  is  fairly 
localized,  serving  to  increase  Kt  for  crack  points  in  the  vicinity  of  the  hole  bore.  The  free  surface  Kt, 
however,  is  only  slightly  affected  as  crack  separation  approaches  zero.  Table  3  lists  the  hole  bore  and 
free  surface  interaction  factors  for  all  aspect  ratios  considered.  The  decrease  in  K  for  the  double 
cracked  case  may  be  caused  by  elimination  of  through-the-thickness  bending  for  the  symmetric  crack 
configuration. 


SUMMARY  AND  CONCLUSIONS 

Several  conclusions  can  be  made  about  stress  intensity  factors  for  two  symmetric  coalescing  corner 
cracks.  First,  as  tw  o  symmetric  corner  cracks  approach  coalescence  (a/T  =  0.5),  the  greatest  increase  in 
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Fig.  11.  Summary  of  the  effect  of  crack  spacing  t„pla  and  crack  shape  ale  on  stress  intensity  factor 
interaction  at  hole  bore  track  location  (<S  =  90’). 
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Fig.  12.  Summary  of  the  effect  of  crack  spacing  !„pfo  and  crack  shape  ale  on  stress  intensity  factor 
interaction  at  front  surface  location  Ut>  -  0)  for  corner  cracked  hole. 
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Table  3  Dimensionless  interaction  factor,  y.  defined  as  ratio  of  stirss  intensity  factor  for  symmetric  double  crack 
geometry  disided  by  single  corner  crack  result  Interaction  factors  given  at  hole  bote  (<j>  =  90't  and  free  surface 
(4>  =  O')  locations  as  functions  of  dimensionless  separation  along  hole  bore  (l„P/o)  and  crack  shape  ale 


interaction  (actors 

HOLE 

BORE  (« 

.  90*) 

J  FREE  SURFACE  («  *  0*) 

«/c 

1.11 

1.5 

2.0 

3.0 

1.11 

1.5 

2.0 

3.0 

w 

.025 

1.352 

1.312 

1.278 

1.228 

.894 

.931 

.957 

.988 

.051 

1.216 

1.169 

1.136 
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.  B9B 

.932 
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1.084 

1.053 
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.934 

.958 

.987 

.222 
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1.035 

.912 
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.500 
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.986 

.982 

.919 

.948 

.964 

.990 

.857 

.983 

.982 

.987 

.995 

.932 

.961 

.979 

1.009 

1.333 

.952 

.955 

.958 

.963 

.957 

.974 

.999 

1.020 

2.000 

.958 

.964 

.979 

.980 

.969 

.989 

1.033 

1.041 

stress  intensity  factor  occurs  at  the  crack  point  along  the  hole  bore.  The  crack  point  at  the  free  surface 
does  not  experience  a  significant  stress  intensity  factor  increase,  but  shows  a  slight  decrease  in  Kt  as  the 
crack  tips  approach  along  the  hole  bore. 

For  crack  points  along  the  hole  bore,  crack  interaction  (stress  intensity  factor  increase)  does  not 
occur  until  the  crack  separation  distance  is  equal  to  or  less  than  half  of  the  crack  length  a.  The  localized 
effect  of  coalescence  is  found  to  depend  strongly  on  crack  size  ( alT )  and  only  weakly  on  crack  shape 
(u/c).  The  crack  shape  dependence  indicates  that  the  deeper  cracks  (smaller  ale  ratios)  have  a  slightly 
larger  K  (maximum  10.7%  difference)  at  the  hole  bore  location  than  more  shallow  flaws  (larger  ale)  with 
the  same  crack  length  a. 
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Coalescence  of  Multiple  Fatigue  Cracks  at  a  Notch 
R.  Perez*  and  A.F.  P.randt,  Jr.** 

This  paper  describes  experimental  results  from  fatipue  tests  con¬ 
ducted  with  semicircular  notched  polymer  specimens  which  contained  two 
cracks  as  shown  schematically  in  Fia.  1.  The  two  corner  cracks  were 
located  at  opposite  sides  of  the  semicircular  notch,  with  the  left 
crack  defined  by  the  lengths  a-j  and  c, ,  and  the  right  crack  by  dimen¬ 
sions  a  and  c  as  shov/n  in  Fig. lb.  In  the  same  figure,  d  represents 
the  crack  separation. 


a.  Four-Point  Bend  Specimen 


b.  Corner  crack  dimensions  before  coalescence 

Tig.  1  Specimen  geometry  showing  location  of  the  camera  and 
schematic  views  of  the  crack  plane. 

Five  experiments  were  completed  and  represented  cases  of  symmetric, 
non  symmetric,  in  plane  and  out  of  plane  cracks.  In  each  test,  the  pre¬ 
cracked  specimen  was  subjected  to  cyclic  bending  and  subsequent  crack 
growth  photographed  through  the  ends  of  the  transparent  specimen.  As 
the  fatigue  cracks  grew  the  spacing  d  decreased  to  zero  and  the  two 
separate  cracks  coalesced  into  a  single  flaw  with  the  cusped  shape 
shown  in  Fig.  lc.  Projected  images  of  the  crack  photographs  were 
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CRACK  LENGTH  (  INCHES 


enlarged,  measured,  and  scaled  to  actual  size  and  the  data  recorded  in 
graphical  form.  Figs.  2  and  3  illustrate  two  representative  graphs. 


Fig.  l(Cont.)  Specimen  geometry  showing  location  of  the  camera 
and  schematic  views  of  the  crack  plane. 


THOUSANDS  OF  CYCLES  (  N  ) 

Fig.  2  Crack  growth  before  coalescence  as  a  function  of  elapsed 
cycles  for  symmetric  out  of  plane  crack  test  T-5. 
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Fig.  3  Crack  growth  after  encounter  and  coalescence  as  a  function 
of  elapsed  cycles  for  nonsymmetric  out-of-plane  crack 
test  T-4. 

It  was  observed  that  interaction  effects  between  cracks  were  more 
pronounced  at  the  crack  bore  dimension  (a)  than  at  the  crack  surface 
dimension  (c).  Interaction  was  characterized  by  an  increased  rate  of 
growth  of  the  bore  crack  length  as  shown  in  Fig.  4.  This  figure  was 
obtained  by  shifting  the  crack  growth  curve  for  the  small  crack  horizon¬ 
tally  until  both  cracks  had  equal  initial  sizes.  In  a  study  of  stress 
intensity  factors  for  interacting  surface  cracks  Murakami  and  Nemat- 
Nasser  [1]  point  out  that  if  the  distance  between  the  centers  of  two 
semi  ell iptical  cracks  is  greater  than  the  addition  of  the  length  of  the 
smaller  crack  and  half  the  length  of  the  larger,  then  the  cracks  will 
not  influence  each  other.  Although  this  criterion  for  crack  interac¬ 
tion  involves  surface  flaws  in  plates,  it  was  applied  successfully  in 
this  study  of  interacting  comer  cracks  along  a  notch. 

Out  of  plane  cracks  grew  past  each  other  without  linkup,  and  flaw 
growth  along  the  notch  stopped.  Crack  growth  perpendicular  to  the 
notch  and  into  the  specimen  interior  continued  so  that  a  flat  crack 
front  formed  in  this  direction.  The  crack  tips  at  the  hole  bore  even¬ 
tually  curved  into  each  other  and  coalescence  occurred.  In  plane 
cracks  simply  coalesced  to  form  a  through- the-thickness  crack. 

After  the  cracks  grew  past  each  other  in  the  out  of  plane  crack 
tests,  the  cusp  dimension  h  grew  very  fast  initially,  but  later  slowed 
down  when  a  uniform  flaw  developed.  Similar,  but  less  pronounced,  cusp 
behavior  was  observed  in  the  in  plane  crack  test. 


3 


Perez  and  Grandt 


CRACK  LENGTH  (  INCHES  ) 


0.0  5.0  10.0  15.0  20.0  25.0  30-0 


THOUSANDS  OF  CYCLES  (  N  ) 

Fig.  4  Comparison  between  the  crack  growth  curves  of  the  small 
and  large  cracks  of  test  T-2. 
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ABSTRACT 

This  paper  describes  an  analysis  scheme  for  predicting  the  fatigue  life  of 
components  which  contain  multiple  cracks.  Fracture  mechanics  techniques 
are  used  to  predict  the  initial  growth  and  coalescence  of  two  or  more  cracks 
located  along  the  bore  of  a  hole.  The  results  of  experiments  with  multiply 
cracked  specimens  are  described  and  are  compared  with  crack  growth  predic¬ 
tions  obtained  by  the  analysis  method. 
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INTRODUCTION 

Fatigue  cracks  frequently  initiate  simultaneously  near  stress  concentrations 
formed  by  notches.  These  cracks  may  grow  independently  for  a  period,  and 
then,  depending  on  their  initial  spacing,  may  join  into  a  single  dominant 
crack  which  causes  eventual  failure.  The  objective  of  this  paper  is  to 
describe  an  experimental  and  numerical  program  directed  at  predicting  the 
growth  and  coalescence  of  multiple  fatigue  cracks  located  at  holes  in 
plates. 

The  particular  geometries  of  interest  are  shown  in  Fig.  1.  Here  2r  is  the 
diameter  of  an  open  hole  in  a  plate  of  thickness  2t.  Initial  surface  and 
corner  flaws,  which  are  described  by  portions  of  ellipses  with  semi-axis 
dimensions  a  and  c,  are  located  along  the  bore  of  the  hole.  In  the  present 
work,  the  two  separate  cracks  are  assumed  to  be  coplanar,  and  a  remote 
tensile  stress  is  applied  perpendicular  to  the  crack  plane.  Various  combi-: 
"nations  of  surface  and  corner  cracks  are  considered.  The  next  section  •• 
presents  a  1  ife  prediction  scheme  for  calculating  the  time  required  for  the 
4 initial  flaws  to  grow  together,  coalesce  into  a  single  crack,  and  extend 
to  final  failure.  The  following  section  then  describes  experiments  with 
multiply  cracked  specimens  conducted  to  evaluate  the  numerical  analysis. 


NUMERICAL  ANALYSIS 


The  goal  here  is  to  describe  an  approach  for  predicting  the  initial  growth 
and  eventual  coalescence  of  separate  cracks  located  along  the  bore  of  a  hole. 
The  geometries  of  interest  are  summarized  in  Fig.  1.  The  semi-axis  dimension 
along  the  hole  bore  is  given  by  a,  while  the  crack  depth  is  specified  by  c. 
Subscripts  1  and  2  refer  to  cracks  1  and  2  respectively.  Crack  tip  locations 
are  defined  by  points  1  thru  6  as  shown  in  Fig.  1. 

’.'hen  subjected  to  an  applied  stress,  cracks  1  and  2  will  grow.  Let  the  fati¬ 
gue  crack  growth  rates  for  crack  tips  1,  2,  ...,  6  be  dl/dN,  d2/dN,  ...» 
d6/dN.  Note  that,  in  general,  these  crack  growth  rates  will  have  different 
values.  The  crack  growth  increments  for  a  specified  number  of  cycles  aN  of 
loading  are,  however,  related.  Assume,  for  example,  that  the  depth  c2  of 
crack  2  extends  a  small  amount  A5  during  the  interval  aN.  Now,  aN  is  given 
by  Eq.  1. 


aw  -  &5 
4N  - 

dN 


(1) 


The  growth  Ai  of  crack  tip  i  is  given  by 


. -  _  di  _  di  a5 
M  dN  *  aN  dN  *  d5 

•  dN 


(2) 


where  i  takes  the  values  1  to  6  as  appropriate. 

( 

Now,  if  the  cyclic  stress  intensity  factors  AKi  are  known  at  crack  tips  1, 
2,  ...,  6,  the  conventional  fracture  mechanics  approach  may  be  used  to 
compute  the  corresponding  fatigue  crack  growth  rates  (at  the  various  crack 
positions)  from  Eq.  3. 


iiJ  *  F(AKi)  (3) 

Here  di/dN  is  the  crack  growth  rate  at  location  i  and  F  (AKi)  is  the  appro¬ 
priate  fatigue  crack  growth  model  for  the  material  of  interest. 

x 

Stress  intensity  factors  for  the  multi-cracked  holes  were  obtained  in  the 
following  manner.  First,  the  Newman  and  Raju  (1981)  solutions  for  single 
surface  or  corner  cracks  located  at  holes  were  used  to  compute  K  at  crack 
tip  locations  1  thru  6.  Newman  and  Raju  (1981)  have  fit  empirical  equations 
to  earlier  finite  element  results  obtained  for  the  single  crack  geometries. 
Their  empirical  solutions  are  readily  coded  for  computer  use,  and  provide 
an  effective  means  for  interpolating  between  the  original  finite  element 
results  obtained  for  discrete  crack  shapes.  - 

The  effect  of  the  adjacent  flaw  was  accounted  for  by  an  "interaction  factor* 
obtained  previously  by  Heath  and  Grandt  (1984).  In  that  work,  the 
three-dimensional  finite  element-  alternating  method  was  used  to  compute 
stress  intensity  factors  for  symetric  corner  cracks  located  on  opposite 
sides  the  plate  thickness  (the  symmetric  version  of  Fig.  lc).  The  double 
crack  K  was  then  divided  by  the  corresponding  result  for  a  single  flaw  to 
obtain  the  dimensionless  interaction  factor  y  given  in  Fig.  2.  Here  the 


INTERACTION  FACTOR 


increase  in  K  at  the  point  where  the  corner  crack  intersects  the  hole  bore 
(locations  3  and  4  in  Fig.  lc)  is  given  as  a  function  of  dimensionless  crack 
spacing  ts/a  and  crack  shape  a/c.  Note  in  Fig.  2  that  as  the  cracks  approach 
coalescence  (ts/a  <0.3),  K  is  increased  signif icantly  along  the  hole  bore. 

In  addition,  flaw  shape  has  a  relatively  minor  effect  on  the  increase  in  K. 
Since,  Heath  and  Grandt  shoved  that  the  stress  intensity  factor  at  free  sur¬ 
face  locations  2  and  5  is  relatively  unaffected  by  the  approaching  crack, 
no  interaction  was  considered  at  those  points,  and  the  single  crack  Ncwman- 
Raju  (1981)  solutions  were  used  to  compute  K  there. 

Thus,  Equations  1  through  3  were  solved  in  an  iterative  manner  which  allows 
the  cracks  to  grow  naturally.  When  the  crack  tips  3  and  4  touch,  a  single 
corner,  surface,  or  through-the  thickness  flaw  was  assumed  as  appropriate, 
and  the  corresponding  single  crack  K  solution  used  to  continue  the  analysis. 
Other  possible  comb i nations,  such  as  penetration  of  surface  crack  locations 
1  or  6  through  a  free  surface  to  form  corner  flaws  were  also  considered  in 
the  algorithm.  Additional  details  are  given  by  Tritsch  (1983).  Sample 
calculations  obtained  by  this  multi-degree-of-freedom  analysis  are  compared 
with  exper iiuontal  results  in  the  following  section. 


LXP1  R1I1LN1 AL  PROCEDURE 

This  section  describes  a  set  of  experiments  conducted  with  multi-cracked 
holes  to  provide  initial  evaluation  of  the  predictive  scheme.  The  double 
corner  crack  geometry  given  in  Fig.  lc  was  considered  here. 

Test  specimens  were  machined  from  a  single  sheet  of  polymethylmethacrylate 
(PMMA),  a  transparent  polymer.  The  specimens  were  25.4  cm  (10.0  in.)  long, 
8.9  cm  (3.5  in.)  wide,  1.9  cm  (0.75  in.)  thick,  and  contained  a  central 
0.95  cm  (0.375  in.)  diameter  hole.  The  specimens  were  annealed  to  relieve 
possible  residual  stresses,  and  the  ends  were  polished  to  transparency. 

Small  notches  were  cut  with  a  razor  blade  along  the  hole  bore  (on  opposite 
sides  of  the  plate  thickness)  to  provide  crack  initiation  sites.  The  speci¬ 
mens  were  then  precracked  in  cyclic  bending  until  fatigue  cracks  had 
completely  surrounded  the  pre-notches. 

Following  precracking,  grips  were  bonded  to  the  ends  of  the  specimens  and 
a  constant  amplitude  tensile  load  was  applied  at  2  Hz  with  a  closed  loop 
electrohydraul ic  fatigue  machine.  Strain  gages  were  mounted  to  the  front 
and  back  sides  of  the  specimens  to  verify  that  uniform  tension  was  applied 
without  bending.  A  cutout  in  one  end. of  the  grips  allowed  a  mirror  to  be 
placed  at  an  angle  over  the  transparent  specimen  end  and  provided  optical 
access  to  the  crack  plane.  Crack  growth  was  recorded  with  a  35  mm  camera, 
and  the  film  strips  measured  to  give  crack  dimensions  as  a  function  of 
applied  load  cycles.  Additional  experimental  details  are  given  by 
Perez  (1933). 

Fatigue  crack  growth  curves  for  two  specimens  are  given  in  Figs.  3  and  4. 
Here  the  experimental  measurements  are  given  by  the  data  points  defined  in 
the  legends  on  Figs.  3  and  4.  Measurements  are  given  for  the  four  crack 
dimensions  al  and  cl  (crack  1)  and  a2  and  c2  (crack  2)  prior  to  coalescence. 
Following  link  up,  the  surface  dimensions  cl  and  c2  of  the  through-the- 
thickness  flaw  are  also  given. 

The  solid  and  dashed  lines  in  Figs.  3  and  4  are  predictions  obtained  by  the 
numerical  scheme  discussed  in  the  last  section.  Baseline  fatigue  crack 
growth  data  (da/dN  versus  aK)  were  obtained  with  additional  through-cracked 


specimens  obtained  from  the  same  sheet  of  FMMA  and  loaded  at  the  same  cyclic 
frequency.  A  simple  power  law  was  used  to  fit  these  baseline  data  and 
provided  the  crack  growth  model  used  for  Fq.  3.  Care  wa s  taken  not  to 
extrapolate  the  fatigue  crack  growth  law  for  the  baseline  data  beyond  its 
1 imits  of  val idity . 

Note  that  the  predictions  of  Figs.  3  and  4  generally  agree  quite  well  with 
the  experimental  data.  In  both  cases,  coalescence  occurred  slightly  before 
the  prediction,  but  actual  fracture  occurred  shortly  following  the  analysis. 
The  latter  result  is  more  likely  due  to  conservatism  in  the  analysis  scheme 
resulting  from  the  fact  that  transition  to  a  uniform  through-the-thickness 
flaw  was  assumed  to  occur  immediately  following  coalescence  (joining  of 
points  3  and  4).  In  actuality,  a  short  period  of  time  was  required  for  the 
joined  cracks  to  grow  into  a  uniform  front.  More  detailed  measurements  of 
the  transition  crack  growth  are  given  by  Perez  (1983). 


C0NCIUD1NG  REMARKS 

The  objective  of  this  paper  has  been  to  describe  a  procedure  for  predicting 
the  lives  of  components  which  contain  two  or  more  adjacent  cracks.  A  multi¬ 
degree  of  f ref  don  model  has  been  described  for  computing  the  initial  growth, 
coalescence,  and  final  extension  of  corner  and  surface  cracks  located  along 
the  Lore  of  a  hole  in  an  open  plate.  Experiments  described  here  with  corner 
cracked  holes  in  transparent  polymer  specimens  gave  excellent  agreement  with 
the  predictive  model.  Similar  predictions  have  been  made  by  the  authors  for 
experiments  in  more  conventional  structural  metals,  and  also  show  good 
agreement  between  experiment  and  analysis.  (It  is  planned  to  describe  those 
results  in  a  separate  paper.) 

Although/ the  present  model  is  limited  to  two  initial  cracks,  the  procedure 
is  readily  extended  to  additional  flaws.  The  stress  intensity  factor  calcu¬ 
lations  for  interacting  cracks,  based  on  the  single  crack  solutions  given 
by  Newman  and  Raju  (1981)  and  modified  by  the  crack  interaction  factor 
obtained  by  Heath  and  Grandt  for  symmetric  corner  cracks,  worked  well  for 
the  experiments  which  have  been  examined  to  date.  It  may  be  necessary  to 
develop  alternate  interaction  factors  for  crack  combinations  which  differ 
greatly  in  size,  as  the  Heath  and  Grandt  result  assumes  symmetric  flaws. 

In  addition,  cracks  which  initiate  in  different  planes  may  need  to  be 
analyzed  by  an  alternate  procedure. 
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ABSTRACT 

i 

The  objective  of  this  paper  is  to  describe  results  of  an 
experimental  and  numerical  study  concerned  with  predicting  the 
growth  and  coalescence  of  multiple  fatigue  cracks.  Naturally 
occuring  fatigue  cracks  frequently  initiate  at  multiple  sites 
near  stress  concentrations.  Individual  cracks  then  grow  by 
cyclic  loading,  until  eventually  two  or  more  adjacent  cracks  link 
up  into  a  single,  dominant  flaw  which  grows  to  final  failure. 

Figure  1  shows  schematic  examples  of  various  surface  and  corner 
crack  combinations  located  at  the  bore  of  a  notch.  The  current 
paper  describes  a  life  prediction  scheme  for  computing  the  total 
life  of  such  multiply  cracked  notches. 

The  fatigue  crack  growth  algorithm  is  based  on  conventional 
linear  elastic  fracture  mechanics  concepts.  Stress  intensity 
factors  are  computed  at  various  points  along  the  crack  tip  (such 
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as  points  1  to  6  in  Fig.  la).  The  fatigue  crack  growth  law  for 
the  component  material  is  then  used  to  compute  an  increment  of 
crack  growth  at  these  flaw  locations.  The  crack  shape  then 
changes,  and  the  procedure  is  repeated  in  an  iterative  fashion 
until  the  flaws  coalesce,  penetrate  a  free  surface,  or  cause 
fracture.  The  cycles  required  for  growth  of  the  multiple  cracks, 
and  the  resulting  single  coalesced  flaw,  are  summed  as  the  itera¬ 
tion  proceeds  to  give  total  life. 

The  stress  intensity  factors  are  based  on  solutions  by  New¬ 
man  and  Raju  (1)  for  single  surface  and  corner  cracks  located  at 
holes.  These  single  crack  results  are  then  modified  by  an 
"interaction  factor"  obtained  from  a  stress  intensity  factor 
solution  for  multiple  corner  cracks  (2).  In  this  manner,  the 
influence  of  the  adjacent  flaw  on  the  stress  intensity  factor  at 
crack  tip  locations  3  and  4  is  determined  (see  Fig.  la).  Once 
crack  tips  3  and  4  join,  a  single  surface  or  corner  crack  is 
assumed,  and  the  appropriate  K  solution  for  the  single  flaw  (1) 
is  employed.  If  the  crack  continues  to  grow,  and  becomes  a 
through-the-thickness  flaw,  then  the  appropriate  through-crack  K 
solution  is  used.  Additional  modifications  to  consider  semicir¬ 
cular  edge  notches  are  described  in  the  paper. 

A  set  of  experiments  with  multiply  flawed  specimens  were 
conducted  to  evaluate  the  multiple  degree-of-freedom  life  predic¬ 
tion  analysis  outlined  above.  Waspalloy  and  titanium  specimens 
were  machined  to  the  configuration  shown  in  Fig.  2.  Surface 
cracks  were  initiated  at  small  elox  starter  slots  located  along 
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the  bore  of  one  of  the  semicircular  edge  notches.  The  specimens 
were  then  cycled  to  failure  under  constant  amplitude  loading 
applied  perpendicular  to  the  crack  plane.  The  fracture  surfaces 
were  periodically  marked  by  heat  tinting  procedures  (the  test 
temperature  was  briefly  changed)  so  that  the  crack  dimensions 
could  be  determined  following  final  specimen  fracture. 

Figure  3  compares  experimental  and  numerical  results  for  a 
typical  test.  The  open  symbols  are  measures  of  crack  dimensions 
obtained  by  the  heat  tinting  technique,  while  the  solid  lines 
represent  predictions  by  the  numerical  analysis.  The  results  of 
approximately  one  dozen  tests  are  presented,  and  show  that,  in 
general,  the  analysis  scheme  provides  an  excellent  estimate  for 
the  initial  growth  of  the  multiple  cracks,  their  eventual  coales¬ 
cence  into  a  single  crack,  and  growth  of  the  final  flaw. 
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Figure  1  Schematic  of  multiple  cracks  located  along  bore  of  a  hole 
in  a  large  plate.  Hole  diameter  is  0,  plate  thickness  is  t 
and  remote  tensile  stress  is  applied  in  z  direction. 


Figure  3  Comparison  of  predicted  and  experimental  growth  and  coalescence  of 
two  surface  cracks  at  semicircular  notch  in  Waspaloy  specimen 


